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1.  Introduction 

Controllability  and  observability  are  two  of  the  basic  questions 
in  control  theory  [5,1 1  ].  Formally,  controllability  denotes  the  ability 
of  steering  a  system  from  a  given  initial  state  to  a  desirable  final 
state;  whereas  observability  is  the  ability  to  determine  uniquely 
the  state  of  the  system  from  observable  quantities. 

The  controllability  of  viscoelastic  fluids  under  shear  flow  was 
studied  in  [12-14]  and  the  controllability  of  viscoelastic  fluids 
under  other  flow  fields  was  investigated  in  [15,16].  Relatively  little 
research  has  been  conducted  on  the  observability  of  viscoelastic  flu¬ 
ids.  However,  observability  is  a  very  important  concept  in  control 
theory  since  it  measures  how  well  internal  states  of  a  dynamical 
system  can  be  inferred  by  knowledge  of  its  external  outputs.  In 
practice,  noise  may  appear  in  both  the  observations  and  the  sys- 
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tern  dynamics.  If  the  system  is  not  observable,  then  we  will  not 
be  able  to  get  an  accurate  estimate  of  the  state,  which  will  defi¬ 
nitely  reduce  (if  not  completely  destroy)  our  ability  to  guide  the 
system  to  a  desirable  final  state.  Therefore,  it  is  critical  to  study  the 
observability  of  a  dynamical  system. 

In  this  paper  we  extend  our  previous  works  to  address  the 
observability  of  viscoelastic  fluids.  We  are  following  earlier  work  on 
Eulerian  and  Lagrangian  observability  of  point  vortex  flows  found 
in  Krener  [8[.  In  all  the  models  considered  here,  we  impose  the 
velocity  field  and  assume  the  viscoelastic  stress  is  homogeneous. 
We  assume  the  first  normal  stress  difference  is  measured  in  exper¬ 
iments.  Our  goal  is  to  infer  the  full  viscoelastic  stress.  The  first 
normal  stress  difference  is  an  important  rheological  property  of 
complex  fluids.  Common  experimental  devices  of  measuring  the 
first  normal  stress  difference  include  cone-and-plate  as  well  as 
parallel-plate  rheometers  [3,6,10].  It  is  also  possible  to  obtain  val¬ 
ues  of  the  first  normal  stress  difference  using  hole  and  exit  pressure 
data  [1], 

We  apply  the  observability  rank  condition  to  study  short-time 
local  observability  of  various  models  from  the  measurement  of  the 
first  normal  stress  difference.  The  mathematical  advantage  of  using 
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the  observability  rank  condition  is  that  it  is  simply  algebraic  and 
systematic  and  it  also  provides  a  tool  to  study  more  complicated 
systems.  To  the  best  of  our  knowledge,  this  work  is  the  first  attempt 
to  study  the  observability  of  complex  fluids. 

The  outline  of  this  work  is  as  follows.  Section  2  gives  a  brief 
overview  of  the  observability  of  dynamical  systems  and  the  observ¬ 
ability  rank  condition.  Sections  3-6  are  devoted  to  the  analysis 
on  the  observability  of  various  viscoelastic  models  under  shear 
or  extensional  flow  by  applying  the  observability  rank  condition. 
Section  7  presents  some  numerical  experiments  where  the  per¬ 
formance  of  the  unscented  Kalman  filter  is  examined  both  in  the 
case  of  zero  observation  noise  and  in  the  case  of  Gaussian  obser¬ 
vation  noise.  The  main  results  of  our  study  are  summarized  in 
Section  8.  Mathematical  details  of  the  unscented  Kalman  filtering 
are  included  in  Appendix  A. 

2.  Observability  rank  condition  of  dynamical  systems 

For  reader’s  convenience,  we  quickly  outline  the  observability 
of  dynamical  systems  and  observability  rank  condition.  In  the  sub¬ 
sequent  sections  we  will  investigate  the  observability  of  several 
viscoelastic  models  using  the  mathematical  tools  introduced  in  this 
section. 

Consider  the  following  system  without  a  control  input 


x=/(x),  xeRn 

(1) 

y  =  h(x),  yeRp 

(2) 

x(0)  =  x°, 

(3) 

where  x  is  the  state  variable,  y  is  the  quantity  that  can  be  observed 
experimentally  and  relation  (2)  is  called  an  observation  equation. 
We  shall  assume/  and  h  are  sufficiently  smooth  functions.  The  state 
x  is  not  observed  directly  but  the  output  y  is.  Can  we  determine  the 
initial  state  x°  from  the  output  history  y(0  :  oo)?  In  this  definition, 
thesymboly(0  :  T)  represents  the  trajectory  t  ->•  y(t)  where  0  <  t  < 
T.  Is  the  map  x°  ->■  y(0  :  oo)  one  to  one?  If  so,  then  mathematically 
we  can  reconstruct  x°  fromy(0  :  oo)and  we  say  the  system  (1 )— (2) 
is  o  bservable.  Note  that  x(t)  can  be  viewed  as  the  initial  state  for 
evolution  beyond  time  t.  Therefore,  the  observability  also  implies 
thatx(t)  can  be  calculated  fromy(0  :  oo). 

There  are  other  possible  definitions  of  observability.  For  exam¬ 
ple,  the  system  (1 )— (2)  is  s  hort-time  observable  if  for  every  T  >  0 
the  map  x°  ->■  y(0  :  T)  is  one  to  one.  In  other  words,  the  effect  of  the 
initial  data  is  immediately  felt  by  the  system.  The  system  is  short- 
time  locally  observable  if  for  every  T  >  0,  the  map  x°  ->■  y(0  :  T)  is 
locally  one  to  one.  Locally  one  to  one  means  that  around  each  x° 
there  is  a  neighborhood  such  that  the  map  is  one  to  one  on  the 
neighborhood.  More  specially,  locally  one  to  one  means  that  given 
any  x°,  there  is  a  neighborhood  U(x°)  such  that  if  x1  e  U(x°),  then 
output  from  x1  is  different  from  that  from  x°.  short-time  locally 
observable  is  the  most  useful  and  easiest  to  measure.  How  do  we 
decide  whether  a  system  is  short-time  locally  observable?  In  the 
following  we  review  a  sufficient  condition  for  the  short-time  local 
observability  [8]. 

First,  we  need  to  introduce  some  notations.  Note  that  from  (2) 
and  (3)  we  havey(O)  =  h(x°)  and  thus 

f(0).g(*W)  (4) 

using  the  chain  rule  and  (1).  The  right  hand  side  of  (4)  is  called  the 
Lie  derivative  of  the  function  h  by  the  vector  field/: 

Lf(h)(x)  =  ||(x]f(x).  (5) 


Notice  that  the  Lie  derivative  is  another  function  from  Rn  to  Rp,  so 
we  can  repeat  the  process  of  taking  the  Lie  derivative: 

Ljmx)  =  Lf(Lf(h))(x)  =  L/(^(x)/(x)) 

L*(h)(x)  =  Lf(Lf-\h))(x). 

Definition  2.1.  gi(x), . . . ,  g/((x)  separate  points  if  given  any  pair 
of  two  points  x°  and  x1,  there  is  at  least  one  g,(x)  such  that 
£i(*°)  A  &(*’)■ 

If  gi(x) .  gkW  separate  points,  then  x -> 

giM\ 

:  is  one  to  one.  Mathematically,  if  the  functions 

gk(x )  / 

ij.  (h)(x),  j  =  0, . . . ,  k,  locally  separate  (or  distinguish)  points 
in  Rn  for  some  k,  then  the  system  is  short-time  locally  observable. 
A  sufficient  condition  for  this  to  happen  is  that  the  one  forms 

dhj(x), . . . ,  dL}  (hj)(x),  i=l,...,p 

span  n  dimensions  at  every  x  where 


Definition  2.2.  The  system  (1 )— (2)  satisfies  the  observabil¬ 
ity  rank  condition  at  x°  if  there  exists  a  k  such  that 

dLj  (h,)  :j  =  0, . . . ,  k;  i  =  1, . . .  ,p|  has  rank  n.  The  system  (1)- 

(2)  satisfies  the  observability  rank  condition  if  it  satisfies  it  at  every 
xeR"  (note  k  may  vary  with  x). 

For  linear  systems,  observability  rank  condition  (ORC)  implies 
global  short-time  observability.  For  nonlinear  systems,  ORC  is  a 
sufficient  condition  of  short-time  local  observability.  Furthermore, 
ORC  is  almost  necessary  for  short-time  local  observability.  Her¬ 
mann  and  Krener  [4]  have  proved  that  if  the  ORC  is  violated  on  an 
open  subset  ofRn,  then  the  system  (l)-(2)  is  not  short-time  locally 
observable. 

3.  The  upper  converted  Maxwell  (UCM)  model 

In  this  section  we  start  with  simple  linear  models.  Then  we 
extend  the  analysis  to  nonlinear  cases  in  the  following  sections. 

A  simple  theory  for  viscoelasticity  was  proposed  by  Maxwell  in 
1 867  [2,9]  and  has  been  usually  called  the  upper  convected  Maxwell 
model.  For  all  the  models  considered  in  this  paper,  we  posit  the 
velocity  field  and  consider  only  the  equation  for  the  stress  tensor. 
In  this  setting,  the  UCM  model  is  linear  and  the  UCM  constitute 
equation  can  be  written  in  the  form 

f-(Vv)T-T(Vv)r  +  AT  =  2/rD,  (7) 

where  T  is  the  stress  tensor,  v  is  the  velocity,  Vv  is  the  velocity  gra¬ 
dient  tensor,  X  is  the  relaxation  rate,  fi  is  the  elastic  modulus  and 
D  is  the  rate-of-deformation  tensor  (i.e.  the  symmetric  part  of  the 
velocity  gradient).  In  the  more  complicated  case  where  the  velocity 
field  is  affected  by  the  stress  tensor,  we  have  a  joint  evolutionary 
system  for  the  velocity  field  and  the  stress  tensor,  which  is  nonlin¬ 
ear.  In  this  paper  we  will  focus  on  the  linear  case  where  the  velocity 
field  is  imposed. 

For  simplicity  we  restrict  our  attention  to  two-dimensional 
homogeneous  viscoelastic  fluids.  We  denote  the  stress  tensor  by 

T-n  T12 

Tn  T22 


(8) 
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The  state  of  the  system  (7)  is  characterized  by  viscoelastic  stress  T 
with  three  components  Tn ,  T22  and  T12. 

Now  we  investigate  the  short-time  local  observability  of  com¬ 
plex  fluids  governed  by  the  Maxwell  equation  under  two  different 
simple  flow  conditions:  extensional  flow  and  shear  flow. 


3.2.  UCM  under  shear  flow 

The  velocity  field  of  a  shear  flow  with  rate  Pe  is  described  by 
v  =  Pe(y,  0)  (17) 


3.1.  UCM  under  extensional  flow 

The  velocity  field  of  a  homogeneous  extensional  flow  with  rate 
Pe  can  be  written  as 

v=y(x,-y),  (9) 

so  the  velocity  gradient  and  the  rate-of-deformation  tensor  both 
have  a  diagonal  form: 


The  rate-of-strain  tensor  is 

0  1 
1  0 

The  UCM  model  (7)  becomes 
tn  —  2PeT12  +ATn  =  0, 
ti2  -  Pe  T22  +  X  T12  =  ij.  Pe, 

t2  2  +  ^t22  =  0. 


1  t  Pe 

D  =  2  [Vv  +  (Vv)T]  =  — 


(18) 


(19) 


„  Pe 
W=y 


1  0 
0  -1 


1  t  Pe 

D=-[Vv  +  (Vv)t]=t 


0 

-1 


(10) 


(11) 


With  these  simplifications  (7)  can  be  expressed  as 
tn  -  (Pe  —  A,)Tii  =  ptPe, 
tl2  +  AT12  =  0, 
t22  +  (Pe  +  X)T2  2  =  —(i  Pe. 

Note  that  the  evolution  of  Tn  and  T22  is  independent  of  Ti2.  So  we 
study  the  subsystem  governing  the  evolution  of  Tn  and  T22: 


We  remark  that  although  the  dynamics  of  T22  can  be  separated 
from  Tn  and  T22,  P\2  still  affects  the  evolution  of  Tn  and  T22:  and 
consequently  Ti2  may  still  be  recovered  from  observations  of  Tn  — 
T22,  which  is  the  case  we  will  see  below. 

The  system  (19)  can  be  rewritten  as 


dx 


dt 
where 


=  /(*), 


(12) 

■*r 

"Til 

’-A.xi+2  Pex2 

X  = 

^2 

= 

Tl2 

,/(x)  = 

-Xx2  +  Pex  3  +  fiPe 

_*3. 

_T2  2. 

-\x3 

(20) 


(21) 


Assume  that  the  quantity  that  can  be  measured  is  the  normal  stress 


i'll  -(Pe-X)TU  =  ix Pe, 
t22  +  (Pe  +  X)r22  =  — p.Pe. 

It  is  convenient  to  introduce  the  following  vector 


(13) 


7(1 

Tn 

7(2 

T22 

This  enables  us  to  rewrite  the  system  (13)  as 


dx  j  _ 

dF=/(x)  = 


7i‘ 

-A.X!  +Pe(/x  +Xi) 

h 

-Xx2  -  Pe(/x  +x2) 

(14) 


Since  the  first  normal  stress  difference  N3  =  Tn  -  T22  can  be  mea¬ 
sured  experimentally  [9],  we  assume  that  the  observation  is  the 
difference  of  Tn  and  T22,  h  =  Tu  —  T22  =  x3  —  x2.  The  Lie  bracket  is 


"3  h  3  h' 

7“ 

3xj  ’  3x2 

h 

3  h. 

~  9xj 

=  [1.-1] 

=  /i  — f2  =  2 )iPe  +  (Pe  -  AJx!  +  (Pe  +  X)x3. 


(15) 


Consequently  we  have 


h(x)  =  Tu  -  T22  =  x-i  -x3. 


To  check  the  observability  rank  condition,  we  first  calculate  the  Lie 
brackets: 

Lj  h  =  dh^(x)  =[1,0,  -1  ]/(x)  =  X(x3  -  )  +  2 Pex2, 

Ljh  =  Lf(Ljh)=[-X,2Pe,X]f(x)  (22) 

=  X2x\  -  4XPex2  +  (2 Pe2  -  X2)x3  +  2/xPe2. 


-  h  - 

'  1 

0  -1 

Lj  h 

= 

-X 

2  Pe  X 

l2  h 

X2 

-4 XPe  2Pe2  -  X2 

L  /  J 

(23) 


This  matrix  is  nonsingular  since  its  determinant  is  4Pe3  =j=.  0.  Thus, 
the  observability  rank  condition  is  satisfied  and  the  Maxwell  model 
(19)  under  shear  flow  is  short-time  locally  observable  if  the  obser¬ 
vation  is  Tn  —  T22.  In  contrast  to  the  case  of  extensional  flow  where 
Ti2  does  not  affect  the  evolution  of  (Tn,  T22)  and  T12  cannot  be 
determined  from  observed  values  of  Tn  -  T22,  in  the  case  of  shear 
flow  all  three  of  Tn,  T22  and  T\2  can  be  reconstructed  from  the 
observed  values  of  Tn  -  T22. 


d 

h(x) 

1  -1 

dx 

Ljh 

Pe  -  X  Pe  +  X 

(16) 


Since  the  determinant  of  this  matrix  is  2 Pe  =/=  0,  the  observability 
rank  condition  is  satisfied.  Therefore,  the  subsystem  of  UCM  model 
under  extensional  flow  is  short-time  locally  observable.  In  other 
words,  when  the  observation  is  Tn  -  T22,  then  Tn  and  T22  can  be 
reconstructed  from  the  observed  values  of  Tn  —  T22.  However,  T32 
cannot  be  determined  from  the  observation  of  Tn  -  T22  since  T12 
does  not  affect  the  evolution  of  Tn  and  T22  at  all. 


4.  The  Phan-Thien-Tanner  (PTT)  model 

The  Phan-Thien-Tanner  model  was  developed  from  network 
theory  to  model  the  rheological  behavior  of  polymer  melts.  It  has 
the  form 

T  -  ( Vv)T  -  T( Vv)T  +  AT  +  v(trT)T  =  2/xD,  (24) 

where  all  the  notations  have  the  same  meaning  as  in  UCM  model. 
In  addition,  the  symbol  “tr”  denotes  the  trace  of  the  tensor  and  v  is 
a  constant. 
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4.1.  PTT  under  extensional flow 


Under  extensional  flow  (9),  the  PTT  model  (24)  can  be  written 
in  vector  form: 


dx 

dt 


=/(*), 


'*1 ' 

'Tn' 

X  = 

*2 

= 

Tn 

.*3. 

_t22_ 

(25) 


and 


"-[A  +  v(xi  +x3)]xi  +(/x  +  Xi  )Pe" 

7i" 

/(*)  = 

— [A  +  v(xi  +x3)]x2 

= 

h 

.-[A  +  v(xi  +x3)]x3  -  {fi  +x3)Pe 

J3  _ 

Notice  that  the  governing  equations  for  x i  and  x3  are  given  by  the 
first  and  third  equations  in  the  above,  which  are  independent  of 
x2.  Therefore,  we  first  consider  the  observability  of  the  subsystem 
consisting  of  two  equations: 


d 

X\ 

_  d 

Tn 

-[A  +  v(xi  +x3)]xi  +(/r  +  Xi)Pe 

dt 

*3 

~  dt 

t22 

-  [ A  +  v(xi  +  x3  )]x3  -  (/x  +  x3  )Pe 

(27) 


2 


Fig.  1.  The  surface  where  the  determinant  in  (31)  is  zero.  The  PTT  model  under 
shear  is  observable  for  all  stress  components  when  the  stress  tensor  is  away  from 
this  surface. 

4.2.  PTT  under  shear  flow 


Again,  we  assume  that  the  observable  quantity  is  the  first  normal 
stress  difference 


In  the  presence  of  a  shear  flow  (17),  the  PTT  model  (24)  becomes 


h  =  Tu  —  T22  =  —  x3. 


Then  it  is  found  that 


det  (J 

h 

|  —  det 

\  OX new 

L-f  h 

.  /new 

J  i 

=  2v(x3  -X! 

+  2Pe, 

1  -1 
2vxi  +  Pe  A  +  2v>X3  +  Pe 

(28) 


where  xnew  = 


.  Clearly,  the  ORC  is  satisfied  ifx3  -  Xj  ^  -  Pe/ v. 


Or  equivalently,  when  the  first  normal  stress  Tn  -  T2 2  does  not 
have  the  value  Pe/v,  then  Tn  and  T22  are  short-time  locally  observ¬ 
able.  That  is  .7,,  and  T22  can  be  reconstructed  from  the  observed 
values  of  Tn  -  T22. 

Now  let  us  briefly  check  the  observability  of  the  whole  system 
(25)  when  the  observation  is  the  first  normal  stress  difference  T11  ~ 
T22.  In  this  case,  the  Lie  bracket  is 


Ljh=  [1, 0, -Ilf  =h-f3. 

Note  that  both/!  and/2  are  independent  of  x2.  It  follows  that  the 
Lie  bracket 


Ljh  =  Lj(Ljh), 

is  also  independent  of  x2 .  As  a  result,  the  second  column  of  the  3x3 
matrix 


d_ 

dx 


Ljh 
L  LIh 


is  entirely  made  of  zeros.  So  the  observability  rank  condition  is  not 
satisfied  and  this  implies  that  from  the  observation  of  Tn  -  T22,  it 
is  not  possible  to  reconstruct  the  whole  viscoelastic  stress  T. 

In  conclusion,  for  the  PTT  model  under  extensional  flow,  the 
observation  of  the  first  normal  stress  difference  Tn  -  T22  allows 
short-time  local  observability  of  Tn  and  T22  when  the  stress  tensor 
is  away  from  Tn  -  T22  =  Pe/v.  T]2  is  not  observable. 


'*1 ' 

'Tn' 

X  = 

*2 

= 

Tn 

.x3. 

T22 . 

and 


'-[A  +  v(Xi  +x3)]xi  +2Pex2 

r/ti 

m= 

—[A  +  v(Xi  +x3)]x2  +  Pex3  +  fjiPe 

= 

h 

.-[A  +  v(xi  +x3)]x3 

h 

If  the  observation  is 


(29) 


(30) 


h  =  Tn  -  T22  =  xi  -X3, 


then  the  Lie  brackets  are 

Ljh=fo-fs  =  -Axi  -  vx2  +  2 Pex2  +  Ax3  +  11X3, 

L?h  =  {-X-  2vx! )/,  +  2 Pef2  +  (A.  +  2vx3  )/3 . 

After  some  algebraic  manipulations,  we  obtain 


=  4Pe  ( Pe 2  -  4vPex2  +  v2x^  -  3v2x|  +  Xvx  1  -  Avx3  +  2v2XiX3)  . 

(31) 


So  one  can  conclude  that  the  PTT  model  under  shear  (29)  is  short- 
time  locally  observable  for  all  stress  components  if  Pe2  -  4  vPex2  + 
v2x2  -  3v2x2  +  Avxi  -  AVX3  +  2u2XiX3  =/=  0. 

In  Fig.  1  we  plot  the  surface  described  by  Pe2  —  4vPex2  +  v2x2  — 
3v2X3  +  Avxi  -  A  vx3  +  2v2XiX3  =  0  where  the  parameters  are  cho¬ 
sen  as  Pe=v  =  A  =  l.  So  the  PTT  model  under  shear  flow  is 
short-time  locally  observable  as  long  as  the  point  (Tn,  T22t  T[2) 
does  not  lie  on  the  surface. 


5.  The  Johnson-Segalman  (JS)  model 

The  Johnson-Segalman  (JS)  model,  which  allows  for  a  non¬ 
monotonic  relationship  between  the  shear  stress  and  shear  rate, 
has  been  widely  used  to  explain  the  striking  “spurt”  phenomenon 
of  non-Newtonian  fluids.  This  phenomenon  describes  a  sudden 
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increase  in  the  volumeric  flow  rate.  The  JS  model  is  given  by 

T-  1  [( Vv)T  +  T(Vv)t]  —  [(Vv)rT  +  T(Vv)]  +  XT=2/xD. 

(32) 

Here  a  is  a  parameter  describing  polymer  slip  where  -1  <  a  <  1; 
When  a  =  1,  the  model  (32)  reduces  to  the  Oldrody-B  model. 

5.2.  JS  under  extensional flow 

Applying  the  extensonal  flow  (9),  we  can  express  the  JS  system 
(32)  as 


This  leads  to 


dx  - 
dt  =S 

and 
/(*)  = 


'*1  ' 

I'll  " 

X  = 

*2 

= 

Tn 

.x3. 

T22 . 

— A.X]  +  {(i  +  ax  i  )Pe 
—Xx2 

—Xx3  —  (/r  +  ax3)Pe 


(33) 


(34) 


As  in  the  upper  convected  Maxwell  model,  the  component  x3  or 
T12  does  not  affect  the  evolution  of  In  and  T22-  So  we  consider  the 
subsystem  of  Tn  and  T22- 


dx  n 


dt 


—  fnew(x)  — 


-Xx-[  +(//.  +  ox  1  )Pe 

Xl 

-A.X3  -  {fi  +  ax3)Pe 

,  Xnew  = 

X3 

(35) 


Assume  the  observation  is  the  first  normal  stress  difference 
h  =  T\\  -  T22  =  x\  - x3, 
then  it  follows  that 


d 

h 

1  -1 

dx 

L-f  h 

L  / new  J 

-X  +  a  Pe  X  +  aPe 

(36) 


whose  determinant  is  the  nonzero  value  2 aPe.  This  means  that  the 
observability  rank  condition  is  satisfied.  Therefore,  for  the  JS  model 
under  extensional  flow  the  observation  of  the  first  normal  stress 
difference  In  -  T22  allows  short-time  local  observability  of  Tn  and 
T22,  but  not  T12. 

5.2.  JS  under  shear  flow 

Under  shear  flow  the  JS  model  is  described  by 


dx 

dt 

and 


/(*)  = 


■*1  ■ 

'Tn' 

X  = 

x2 

= 

Tn 

_x3  _ 

T22 . 

(37) 


-Xxi  +  (a  +  l)Pex2 

,  a  +  1  n  a  -  1  „ 

-Xx2  h - 2 — Pex3  H - 2 — Pexi  +  M Pe 


'ft' 

f2 

h 

(38) 


—Xx3  +  (a  -  l)Pex2 
If  the  observation  is 
h  =  fii  ~  T22  =  - x3, 

then  the  Lie  brackets  are 
Ljh  =fi  —f3  =  -Xxi  +  Xx3  +  2  Pex2, 

Lj  h  =  -A/,  +  X/3  +  2 Pef2 

=  [X2  +  (a  -  l)Pe2]xi  —  4XPex2  +  [-X2  +  (a  +  l)Pe2]x3  +  2/xPe2. 


det|  95 


Ljh 


=  det 


Llhd 

1  0  -1 

—X  2  Pe  X 

X2  +  (a  -  1  )Pe2  -4 XPe  -X2  +  (a  +  1  )Pe2 


=  4aPe3  =f=  0. 

(39) 


So  the  observability  rank  condition  is  validated  and  therefore  for 
the  JS  model  under  shear  (37)  the  viscoelastic  stress  is  short-time 
locally  observable  if  the  observation  is  the  first  normal  stress  dif¬ 
ference. 

6.  The  Giesekus  model 

The  Giesekus  model  is  a  nonlinear  viscoelastic  fluid  model  which 
reproduces  many  characteristics  of  the  rheology  of  polymer  solu¬ 
tions  as  well  as  other  liquids.  It  takes  the  form 

T  -  (Vv)T  —  T(Vv)r  +  >.T  +  r>T2  =  2/xD.  (40) 

6.2.  The  Giesekus  model  under  extensional  flow 

Under  extensional  flow  (9),  the  Giesekus  model  (40)  becomes 
dx 


dt 

and 

/(*)  = 


'xt ' 

'Tu' 

X  = 

x 2 

= 

Tn 

_x3  _ 

T22 . 

'  -A.X!  -  v(x2  +  x2)  +  (/r  +  X!  )Pe 

'ft' 

-[X  +  u(x,  +x3)]x2 

= 

f2 

— XX3  -  v{xfj  +  x| )  -  (/r  +  X3  )Pe 

h 

(41) 


(42) 


To  investigate  the  observability  of  the  system  (41 ),  we  assume  the 
observation  is  the  first  normal  stress  difference 

h  =  T11  —  I22  =  x,  -  X3. 

Then,  after  some  calculations,  we  find  that 

h 


Ljh 

L2h 

■  f 


=  4vx2[v(x1  -  x3)  —  Pe]2 


(43) 


So  the  observability  rank  condition  test  is  valid  if  x2  0  (i.e. 
T12  =/=  0)  and  v(x3  -  x^ )  +  Pe  0  (i.e.  Tn  -  T22  Pe/v).  That  is,  for 
the  Giesekus  model  under  extensional  flow,  the  observation  of  the 
first  normal  stress  difference  Tn  —  T22  gives  the  short-time  local 
observability  of  the  viscoelastic  stress  provided  that  T12  i=  0  and 
I’ll  -  T22  +  Pe/v. 

6.2.  The  Giesekus  model  under  shear  flow 

Finally,  we  consider  the  observability  of  the  Giesekus  model 
under  shear  flow: 


dx 


dt 
and 

/(*)  = 


=  /(x),  x  = 


'xt ' 

'Tu' 

X2 

= 

Tn 

.x3. 

T22 . 

— Xx^  -  v(x2  +  x%)  +  2 Pex2 
-[A  +  v(xi  +x3)]x2  +  Pex3  +  fiPe 
—Xx3  —  v(x^  +  x ^ ) 


(44) 


(45) 
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Fig.  2.  The  surface  where  the  determinant  in  (46)  vanishes.  The  Giesekus  model 
under  shear  is  observable  for  all  stress  components  as  long  as  the  stress  tensor  is 
away  from  this  surface. 


As  usual,  assume  that  the  observation  is  the  first  normal  stress 
difference 

h  =  Tn  -  T22  =  *i  —  *3- 
Some  algebra  gives 

^  =  4 Pe3  +4PeXv(xi  -X3)-  16v(v2x1X3  +Pe2)x 2 
+8v2Pe(xi  -  X3  )x3  +  8v3X2(x2  +x|). 

(46) 

Therefore,  the  Giesekus  model  under  shear  (44)  is  short-time 
locally  observable  if  the  determinant  in  (46)  does  not  vanish. 

Fig.  2  depicts  the  surface  where  the  determinant  in  (46)  is 
zero.  Here  we  choose  all  the  parameters  to  be  one:  Pe  =  X  =  v  =  1. 
Clearly,  the  Giesekus  model  under  shear  flow  is  short-time  locally 
observable  provided  that  (Tn,  T22,  T12)  is  not  on  this  surface. 

7.  Unscented  Kalman  filtering  (UKF)  of  complex  fluids 

The  observability  of  a  system  tells  us  only  the  mathematical 
possibility  of  recovering  the  state  of  the  system  from  observable 
quantity.  It  does  not,  however,  provide  us  a  concrete  way  of  extract¬ 


Ljh 

T?  h 
L  f  . 


ing  the  state  from  observations.  The  most  intuitive  way  of  direct 
reconstruction  of  state  is  to  use  repeated  differentiation.  Unfor¬ 
tunately,  numerical  differentiation  is  susceptible  to  measurement 
noise.  The  family  of  Kalman  filters  is  an  efficient  and  robust  method 
for  reconstructing  the  solution  without  resorting  to  numerical  dif¬ 
ferentiation.  Numerical  differentiation  requires  only  data  over  very 
short  time.  In  contrast,  Kalman  filters  use  data  over  a  period  of  time 
to  avoid  the  instability  associated  with  repeated  numerical  differ¬ 
entiation.  In  addition,  Kalman  filters  can  deal  with  the  case  where 
the  measurements  are  polluted  by  noise  of  fairly  large  magnitude, 
in  which  case  numerical  differentiation  will  certainly  fail. 

In  this  section  we  demonstrate  the  performance  of  unscented 
Kalman  filtering  in  recovering  the  viscoelastic  stress  from  the 
observed  first  normal  stress  difference.  We  will  examine  two  cases: 
(1)  the  observation  is  noise  free  and  (2)  the  observation  contains 
Gaussian  noise. 

The  unscented  Kalman  filtering  was  first  developed  by  Julier  et 
al.  [7]  and  it  has  been  one  of  the  workhorses  of  nonlinear  estima¬ 
tion  problems.  While  the  original  Kalman  filter  was  developed  for 
linear  systems,  the  idea  was  generalized  to  nonlinear  systems  and 
led  to  many  methods  for  nonlinear  systems,  including  extended 
Kalman  filter  (EKF),  UKF,  Ensemble  Kalman  filter,  etc.  The  technique 
of  UKF  is  more  accurate  and  easier  to  implement  than  an  EKF.  EKF 
is  based  on  the  linearization  of  the  dynamic  system,  which  may 
cause  numerical  instability  due  to  nonlinear  effects.  In  addition, 
EKF  requires  deriving  a  Jacobian  matrix  which  is  a  challenging  task 
for  complex  systems.  In  contrast,  UKF  is  an  approach  that  takes  into 
account  the  nonlinear  dynamics,  rather  than  its  linearization,  in  the 
propagation  of  covariance  matrix.  UKF  is  “founded  on  the  intuition 
that  it  is  easier  to  approximate  a  probability  distribution  than  it  is 
to  approximate  an  arbitrary  nonlinear  function  or  transformation” 
[7].  A  numerical  recipe  for  the  UKF  is  given  in  Appendix  A. 

To  test  UKF  for  a  non-steady  state  solution  of  the  UCM  model 
under  shear  flow,  we  add  an  external  force 


'0 

0.5  sin(t) 
0 


(47) 


to  the  right-hand  side  of  the  equation  (20).  We  remark  that  the 
observability  study  in  the  previous  section  is  still  valid  here  since 
in  (23 )Lj  h  is  unchanged  whereas  3/3 xL?  h  is  the  same  as  before. 

We  start  the  system  (20)  with  the  external  force  (47)  at  some 
initial  condition  x(0)  and  the  filter  at  a  different  initial  condition 
x(0).  We  solve  the  system  equations  without  noise  to  get  state  x(0  : 


Fig.  3.  (a)  Exact  solutions  (solid  lines)  and  estimated  solutions  from  UKF  (symbols)  of  UCM  model  under  shear  flow  with  an  external  force,  (b)  Corresponding  filter  errors  of 
UKF  in  (a). 
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Fig.  4.  (a)  Exact  solutions  (solid  lines)  and  estimated  solutions  from  UKF  (symbols)  of  PIT  model  under  shear  flow  with  an  external  force,  (b)  Corresponding  filter  errors  of 
UKF. 


oo )  and  observation  h( 0  :  oo)  trajectories.  In  the  case  of  observations 
with  no  noise,  we  pass  the  noise  free  observation  trajectory  to  the 
filter  and  the  filter  yields  a  state  estimate  trajectory  x(0  :  oo).  The 
estimation  error  x(t)  -  x(  t)  is  the  difference  between  the  state  of  the 
system  (which  is  not  directly  measurable)  and  the  estimate  state 


produced  by  the  filter  from  the  observation.  The  filter  is  said  to  be 
convergent  if  the  estimation  error  goes  to  zero  as  t  oo  for  any 
x(0)  and  x(0). 

In  Fig.  3  we  show  the  results  of  UKF  for  the  UCM  model  under 
shear  flow.  The  parameters  used  here  are  X  =  1.0,  Pe  =  1.0,  and 


Fig.  5.  (a)  Observation  with  Gaussian  noise  (solid  line)  vs  observation  without  noise  (dashed  line),  (b)  Exact  solutions  (solid  lines)  and  estimated  solutions  from  UKF  (symbols) 
of  UCM  model  under  shear  flow  with  an  external  force  where  observations  contain  Gaussian  noise,  (c)  Corresponding  filter  errors  of  UKF  in  (b). 
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Fig.  6.  (a)  Observation  with  Gaussian  noise  (solid  line)  vs  observation  without  noise  (dashed  line),  (b)  Exact  solutions  (solid  lines)  and  estimated  solutions  from  UKF  (symbols) 
of  PTT  model  under  shear  flow  with  an  external  force  where  observations  contain  Gaussian  noise,  (c)  Corresponding  filter  errors  of  UKF  in  (b). 


yu,  =  1.0.  The  left  panel  depicts  the  exact  and  estimated  solutions 
of  the  UCM  system,  respectively,  whereas  the  right  panel  gives 
the  corresponding  errors  of  the  estimated  solutions.  Fig.  3  shows 
the  convergence  of  the  UKF  when  there  is  no  noise  but  only  initial 
estimate  error. 

Then  we  carry  out  the  simulations  of  the  PTT  model  under  shear 
flow  with  the  external  force  (47)  added  to  the  system  as  well.  Our 
results  are  plotted  in  Fig.  4  where  the  parameters  are  A  =  1 .0,  Pe  = 
1 .0,  /r,  =  1 .0,  and  v  =  1 .0.  Again,  the  UKF  converges  in  the  absence 
of  noise. 

In  reality  all  observations  are  polluted  by  noises  from  various 
sources.  Next,  we  examine  the  performance  of  UKF  in  the  presence 
of  observation  noise.  We  add  Gaussian  noise  to  the  measurements. 
Specifically,  we  assume  that  the  experimental  measurement  is  the 
true  value  of  the  observation  plus  Gaussian  noise: 

/iexperiment(0  :  oo )  =  h(0  :  oo)  +  Gaussian  noise. 

Then  we  pass  the  observation  with  noise  fiexpenment(0  :  oo)  to  the 
filter.  The  filter  outputs  an  estimated  state.  Fig.  5  shows  the  results 
of  UKF  for  the  UCM  model  under  shear  flow.  Fig.  5(a)  depicts  the 
measurement  with  noises  (the  solid  line)  vs  the  measurement  with¬ 
out  noise  (dashed  line).  The  exact  solutions  and  estimated  solutions 
from  UKF  are  plotted  in  Fig.  5(b)  whereas  the  corresponding  filter 
errors  are  given  in  Fig.  5(c).  It  is  clear  that  UKF  gives  very  good 
estimate  even  when  noises  are  present  in  the  measurements. 


The  performance  of  UKF  in  the  presence  of  noise  for  the  PTT 
model  under  shear  flow  is  demonstrated  in  Fig.  6.  Again,  UKF  gives 
very  good  estimates  when  the  measurements  contain  noise.  In  par¬ 
ticular,  the  estimation  error  is  much  smaller  than  the  magnitude  of 
the  measurement  noise. 

8.  Concluding  remarks 

In  this  paper  we  have  applied  the  observability  rank  condition 
(ORC)  to  the  vector  fields  in  the  various  constitutive  models  to 
study  the  short-time  local  observability  of  complex  fluids  driven 
by  extensional  or  shear  flow  fields.  In  each  case  the  measurement 
is  assumed  to  be  the  first  normal  stress  difference  Tn  —  T22.  We 
summarize  our  main  results  as  follows. 

•  For  the  upper  converted  Maxwell  model  under  extensional  flow, 
the  observation  of  the  first  normal  stress  difference  allows  the 
short-time  local  observability  of  the  two  components  Tn  and 
T22  of  the  viscoelastic  stress,  but  not  the  component  T 12.  In  con¬ 
trast,  for  the  UCM  model  under  shear  flow,  the  observation  of 
the  first  normal  stress  difference  leads  to  the  short-time  local 
observability  of  all  components  of  the  viscoelastic  stress. 

•  For  the  Phan-Thien-Tanner  model  under  extensional  flow,  the 
observation  of  the  first  normal  stress  difference  provides  short- 
time  local  observability  of  Tn  and  T2 2,  but  not  T12,  if  Tn  — 
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T22  #  Pe/v.  For  the  PTT  model  under  shear  flow,  the  observation 
of  the  first  normal  stress  allows  the  reconstruction  of  viscoelastic 
stress  as  long  as  the  stress  is  away  from  the  surface  prescribed  by 
the  equation 

Pe2  -  4 vPeTu  +  v2!2,  -  3 v2T\2  +  XvTu  -  XvT22  +  2v2TuT22  =  0. 

•  For  the  Johnson-Segalman  (JS)  model  under  extensional  flow, 
the  observation  of  the  first  normal  stress  difference  leads  to  the 
short-time  local  observability  of  the  two  components  Tu  and 
T22  of  the  viscoelastic  stress,  but  not  the  component  T\2.  For  the 
JS  model  under  shear  flow,  the  observation  of  the  first  normal 
stress  difference  gives  the  short-time  local  observability  of  the 
full  viscoelastic  stress. 

•  For  the  Giesekus  model  under  extensional  flow,  the  observa¬ 
tion  of  the  first  normal  stress  difference  provides  the  short-time 
local  observability  of  the  full  viscoelastic  stress  if  T12  ^  0  and 
Tii  -  T22  =f=  Pe/v.  On  the  other  hand,  the  Giesekus  model  under 
shear  flow  is  short-time  locally  observable  as  long  as  (In,  T22, 
T\2)  satisfies 

4Pe3  +  4PeXv(T„  -  T22)  -  \6v(v2TuT22  +  Pe2)T12 
+  8 v2Pe(Tu  -  T22)T22  +  8v3T12(T121  +  7*)  +  0. 

Finally,  we  point  out  that  all  of  the  four  models  have  observ¬ 
ability  for  all  stress  components  almost  everywhere  under  shear 
flow.  In  contrast,  under  extensional  flow  most  of  the  models  have 
no  observability  for  the  shear  stress  component. 


from  which  the  mean  and  covariance  of  the  estimate  can  be  recov¬ 
ered.  UKF  is  a  filter  where  the  true  mean  and  covariance  can  be 
more  accurately  approximated.  It  can  be  shown  that  the  nonlinear 
transformation  of  the  sigma  points  preserves  statistics  up  to  second 
order  in  a  Taylor  series  expansion.  Based  on  this  fact,  a  prediction 
of  the  state  and  the  covariance  matrices  in  the  filter  algorithm  can 
be  carried  out  as  follows: 

•  From  the  previous  estimate  of  the  state,  xk_] ,  and  the  covariance 
matrix,  P“r  calculate  a  set  of  sigma  points  as 

(%-i,  i  =  l,  2,  ...,n), 

where  (■),■  is  the  i  th  column.  A  sigma  point  is  denoted  by  cr!,  i  = 
0, 1 , . . . ,  2n  and  cr°  =  xf(_1 . 

•  Propagate  all  the  sigma  points  cr,-  through  the  nonlinear  functions 
for  one  step  size  to  obtain  z'.  Update  the  output 

gl  =  Ji(z'),  i  =  1,  2, ...  ,n. 

•  Calculate  the  mean  (prediction)  of  the  state  and  output, 

2  n 

h = y ^w'sz\ 

1=0 
2  n 

yk  = 

i=0 


Acknowledgment 

This  research  was  supported  in  part  by  the  Air  Force  Office  of  Sci¬ 
entific  Research  grant  F1ATA06313G003,  the  Army  Research  Office 
and  the  National  Science  Foundation.  The  authors  thank  the  anony¬ 
mous  referees  for  their  constructive  suggestions  on  improving  this 
manuscript.  The  authors  also  thank  Qi  Gong  for  helpful  discussions 
on  UKF. 

Appendix  A.  The  unscented  Kalman  filtering 

For  systems  with  observability,  Kalman  filter  is  a  widely  used 
tool  to  estimate  the  states  using  observation  data.  In  this  appendix 
we  give  a  short  discussion  on  the  unscented  Kalman  filtering  (UKF) 
whereas  a  detailed  description  of  the  extended  Kalman  filtering 
(EKF)  can  be  found  in  [8], 

Consider  a  system  with  output 

x=/(x),  xeK"  ,  , 

y  =  h(x),  yeRP  v  1 

The  primary  goal  of  Kalman  filter  is  to  estimate  the  unknown  state 
x  using  the  observation  data  ofy.  The  initial  guess  or  estimate  of  x 
is  treated  as  a  random  variable  with  zero  mean  and  a  covariance 
matrix.  It  can  be  shown  that  without  sensor  error,  the  estimate,  x, 
from  Kalman  filter  asymptotically  approaches  the  accurate  value  of 
x.  In  the  presence  of  noise,  the  Kalman  filter  is  able  to  optimize  the 
estimate  process  by  finding  the  estimate  with  smallest  error  (i.e. 
error  with  smallest  standard  deviation). 

Like  the  EKF,  the  UKF  is  derived  by  adding  driving  and  observa¬ 
tion  noises  to  the  above  observed  system.  In  the  UKF,  it  is  assumed 
that  the  estimate  state  x  is  always  a  normally  distributed  variable  at 
every  sampling  instance.  The  mean  and  covariance  information  of 
this  random  variable  can  be  stored  in  a  set  of  specially  selected 
sample  points  called  sigma  points.  These  sigma  points  are  eas¬ 
ily  obtained  from  the  method  described  in  Julier  et  al.  [7].  These 
sigma  points  are  then  propagated  through  the  nonlinear  dynamics, 


where  WJ  is  the  weight  [7]. 

•  The  prediction  of  the  covariance  matrices  is  given  by, 


2  n 

P“  =  y] W‘(z''  -  xk)f  -  xkf, 

i= 0 
2  n 

if  =  y]  wic(,gl -ykXg‘  -faf, 

1=0 
2  n 

if  =  yw^-W-hf, 

i=0 

where  W'c  is  the  weight  for  covariance  matrices  [7].  Here  is 
the  covariance  matrix  of  the  state  and  if  is  the  measurement 
noise  covariance  matrix. 

•  Use  the  predicted  value  of  xk,  P“,  Pf  and  P|*y  to  update  the  esti¬ 
mations  as  follows 


xk=xk  +  Kk(yk  -  yk). 


P?  =  P?-Kkiflf 


where  Kk  is  the  Kalman  gain  and  its  value  is  given  by 
Kk  =  if[iff'.  (49) 

For  a  complete  description  of  the  unscented  Kalman  filtering, 
we  refer  the  readers  to  the  paper  by  Julier  et  al.  [7]. 
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